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ABSTRACT 

One of the most sought-after signatures of reionization is a rapid increase in the ionizing 
background (usually measured through the Lya optical depth toward distant quasars). Con- 
ventional wisdom associates this with the "overlap" phase when ionized bubbles merge, al- 
lowing each source to affect a much larger volume. We argue that this picture fails to describe 
the transition to the post-overlap Universe, where Lyman-limit systems absorb ionizing pho- 
tons over moderate lengthscales ( < 20-100 Mpc). Using an analytic model, we compute the 
probability distribution of the amplitude of the ionizing background throughout reionization, 
including both discrete ionized bubbles and Lyman-limit systems (parameterized by an atten- 
uation length). We show that overlap does not by itself cause a rapid increase in the ionizing 
background or a rapid decrease in the mean Lya transmission toward distant quasars. More 
detailed semi-numeric models support these conclusions. We argue that rapid changes should 
instead be interpreted as evolution in the attenuation length itself, which may or may not be 
directly related to overlap. 
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1 INTRODUCTION 

In the last several years, the cosmological community has made 
an enormous effort to measure the reionization history of the 
intergalactic medium (IGM). But the picture remains murky: 
some evidence (principally the cosmic microwave bac kground, or 
CMB ) points toward a mostly-ionized universe at 2; > 9 jPage et aU 
[2007 : Komatsuetal. 2008; Dunklev et al. 2008l),l)ut other ob- 
servations (principally from quasar absorption) a re usually taken 
to imply that reionization ends only at z ~ 6 jpan et alj 20021: 
Mesinger & HaimanI |2004| : iFan et alj l2006l : iMesinger & HaimanI 



requires a small neutral fraction. Second, there appears to be a sub- 
stantial steepeni ng of the amount of absorption beyond z ~ 6 
( lFanetaljr2006h . However, th e latter conclusion has been chal- 
lenged on empirical grounds jSongaila & Cowid |2002| : ISongailal 
I2OO4I : lOh & Furlanettd [200^ : iBecker et alj l2007h : unfortunately. 



2007h . Many other techniques cannot yet distinguish between 
early and lat e reionizati on s cenarios (e.g., Totani et al. 2006 : 
iKashikawa et a l. 2006; O ta et al.l2007i : lMcOuinn et alj2007il200a : 
see also lpan et al4.2006 for a review). 

While there is not necessarily a contradiction between the 
CMB and quasar measurements - reionization may simply be 
relatively extended, which would hardl y be surprising from a 
theoretical standpoint (e.g.. ICenI I2OO3I: IWvithe & Loebl l2003l : 
iFurlanetto & Loeb|[200^ : llliev et al.ll2007h - the tension does call 
for a critical examination of the existing evidence. The most 
widely-recognized point in favor of late reionization is the appar- 
ent rapi d decrease in the mean Lyct transmission toward quasar s 
at « > 6 jpan et al.ll200lll2002l : lwhite et al.ll2003|:lFan et akbOOfih. 
There are two aspects to this. First, a complete Gunn & Petersonl 
absorption trough appears toward some quasars, although 
the enormous optical depth of that line means that this still only 
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known lines o f sight are sparse enough that no clear resolution has 
emerg ed fe.g., lLidzetalj2006l) . 

Here, we take a complementary approach and examine the the- 
oretical underpinnings of the conclusion that such a rapid change 
implies a detection of "reionization." The crux of this argument 
is that the moment of "overlap" (i.e., the point at which ionized 
bubbles merge into much larger units, usually considered to be the 
"end" of reionization) must be accompanied by a sudden, rapid in- 
crease in the amplitude V of the ionization rate. Qualitatively, this 
expectation comes from percolation models of the reionization pro- 
cess: when ionized bubbles (presumed to be transparent to ionizing 
photons) overlap, many more sources illuminate any given patch, 
so r should increase rapidly. Quantitatively, the first self-consistent 
simul ation of reionization showed precisely such a jump ( Gnedi^ 
I2OOOI) . In that simulation, V evolved much more sluggishly both be- 
fore overlap (when ionized regions grew slowly because the pho- 
tons had to ionize fresh material) and afterwards (when the ionizing 
photons were already able to propagate large distances), so overlap 
appeared to be clearly defined. 

However, this picture does not match properly onto the well- 
understood post-reionization Universe, where dense "Lyman-limit 
systems" (LLSs) absorb ionizing phot ons over relatively small dis- 
tances (e.g., ~ 110 Mpc at 2: = 4; Storr ie-Lombardi et al]|l994 
iMiralda-Esciidi 12003. . probably falling to < 30 Mpc by z = 6; 
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iLidz et alj IiOOtI) . Once ionized regions grew beyond the mean 
separation of tliese systems, LLSs (rather than the edges of ion- 
ized bubbles) controlled the mean free path of ionizing photons 
jpurlanetto & OlJIlOOSljGnedin & Faj|2006h and by extension the 
effective horizon to which sources could be seen. Because reioniza- 
tion is so inhomogeneous, many ionized re gions will reach this size 
( > 20 Mpc) well befo re complete overlap 7Barkana & Loeb" 20041: 
iFurlanetto et al.ll2004[) : there must in fact be a gradual transition 
from the "bubble-dominated" ionization topology characteristic of 
reionization to the "web-do minated" topology cha racteristic of the 
post-reionization Universe jpurlanetto & OhIuOOSi) . 

Here we ask whether overlap must be accompanied by a rapid 
increase in T and, conversely, whether such an increase is a robust 
signal of overlap. We also examine the implications of recent reion- 
ization models for the Lya forest observations conventionally used 
to argue for late reionization. Unfortunately, numerical simulations 
do not yet have the dynamic range to sample the large (~ 100 Mpc) 
scales necessary for reionization and simult aneously predict de- 
tailed properties of the Lya forest (though see lGnedin & Fanll2006l 
for a detailed study of LLSs and the Lya forest during reioniza- 
tion in a 8/i^^ Mpc box). Thus we will use a simple analytic model 
that incorporates both the discrete ionized bubbles and intervening 
absorption by LLSs to show that, during the end stages of reion- 
ization, r is primarily controlled by the mean separation of LLSs, 
which may not evolve rapidly at the moment of overlap - and, more 
importantly, may continue evolving long afterwards. 

In our numerical calculations, we assume a cosmology 
with = 0.26, ST^A = 0.74, Of, = 0.044, Ho = 

100/ikms"^ Mpc"^ (with h = 0.74), n = 0. 95, and gg = 
0.8, consistent with the most recent measurements jPunklev et al.l 
I2OO8. : Komatsu et al. 2008). Unless otherwise specified, we use co- 
moving units for all distances. 



2 METHOD 

We wish to compute the probability distribution /(J) of the angle- 
averaged specific intensity of the radiation background at the HI 
ionization edge, J. Its value at a given point depends on four basic 
parameters: the number density of ionizing sources, n^; the charac- 
teristic luminosity of each source, ; the size of the local ionized 
bubble, R; and the attenuation length for ionizing photons within 
each bubble, ro. This last quantity is determined by the spacing of 
LLSs in the mosfly ionized IGM ( Miralda -Escude 2003). 

We first suppose that R is prescribed and compute / (J) within 
that single ionized region. For the sources, we assume that every 
dark matter halo in the region with a virial temp erature above the 
minim um level for atomic cooling (Tvir > 10* K: lBarkana & Loebl 
1 200 II and references therein) hosts a single galaxy with luminosity 
proportional to its mass, although the mass function is steep enough 
that the luminosity-mass relation makes little difference to our re- 
sults. Our qualitative results are also unaffected if only a fraction 
of galaxies actively form stars (and hence do not produce ionizing 
photons) at any given time; as we will see below, random fluctua- 
tions are small, even if most galaxies are quiescent. 

For our purposes, we can divide the ionized bubbles into two 
limiting regimes. Small bubbles (usually appearing early in reion- 
ization) have i? ^ ro. In this case, we can let ro — + 00 (or, equiv- 
alently, the optical depth t = r /ro experienced by a photon within 
that b ubble vanishes). We then compute f^^o[J) following IZud 
( Il992h . with a simple extensio n to account for the luminosity func- 
tion of the sources (as in iMeiksin & Whit&,2004,) . Assuming that 



the total number of sources in the bubble is Poisson-distributed with 
mean iV — (47r/3)ni_R^_(a reasonable approximation according to 
N-body simulations : ICasas -Miranda et al. 2002), 



/r=0(j) 



ds exp TV j dx (j){x)Ki{sx) 
-sj + N / dx (j>{x)K2{sx) 



(1) 



Here j = J/j;=°, JI=° = L*/(47ri?)2, L„ is the mean lumi- 
nosity of the sources, x = L/ Li,, 4'{x) is the source luminosity 
function normalized so that J dx(j){x) — 1, 



and 



cos(t) - 1 - 2tlmg{t) 
sm{t) +2tRegit) 



du e 



(2) 
(3) 



(4) 



In this limit, the mean background is (j)^^o ~ (J) ^ 

iVJJ^" oc Li,R. 

In the later stages, most bubbles become much larger than the 
attenuation length set by embedded L LSs (ro <^ R), which then 
absorb most of the ionizing photons jFurlanetto &"ohll2005h . In 
this case, we can ap proximate the ionizin g background by taking 
the limit R ^ 00 aMeiksin & Whitell2004h . which gives 



sNo j dx x(f){x)lmG{sx) 



da exp 



sj' + sNo j dx x(f){x)KeG{sx) 



,(5) 



where iVo = (47r/3)n,rg, / = J / J^=°° , Jf=°° = L*/(47rro)^ 



Git) = / diiT^(«)e'*", (6) 
Jo 

and u = e~^/r^. In this limit, the mean background is (j')^^^ — 
3jVo , or (J) oc NoJ^^°° x L^rp. 

iMeiksin & White! ( l2004l) describe how to obtain the full distri- 
bution for arbitrary R and ro; however, we find that the two limiting 
cases above are reasonable approximations to the true distribution 
(and much simpler to compute), provided that they are renormal- 
ized to have the proper mean value that includes both attenuation 
and the finite bubble size jMeiksin & Whitell2004l) : 

(j')=37Vo(l (7) 

The approximation is effective because reasonable attenuation 
lengths at the end of reionization (_R > 10 Mpc) contain huge num- 
bers of sources, for which the limiting forms converge to each other. 
For a given bubble, we therefore use equation l[T) if i? < ro or 
equation l[5) otherwise, rescaled by a constant factor so that f{j) 
has the proper mean. 

Figure \T\ shows some example distributions for bubble s at 
z — 6.5 (with rii computed from the lPress & Schechteij|l974l dis- 
tribution). Note that we have normalized J to its mean value in a 
fully ionized IGM (taking R 00). Clearly the distributions nar- 
row rapidly as more sources become visible (i.e., as R increases); 
this is because the fractional variation in the source counts goes like 
gnoring absorption). However, the peak (and 
the mean) evolve much less rapidly. When R <ti ro, (J) oc R, but 
the evolution slows once 7? ~ ro. Indeed, for a constant ro, the 
mean amplitude increases by only a factor ^ 1.6 from R = ro to 
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V<J>ionLz.d V<J>ion„ed 



Figure 1. Distribution of J relative to its mean value in a fully-ionized 
IGM, (>^)ioniz(,d- curves assume z = 6.5 and tq = 20 Mpc. 
The solid curve is for a fully-ionized universe; the others assume dis- 
crete bubbles with i? = 2, 5, 10, and 20 Mpc, from left to right, or 
N = (0.155, 0.241, 1.93, 15.4) x IQ-*. 

R = oo, as the additional sources are already significantly attenu- 
ated. 

There are three shortcomings of this simple model. First, our 
rescaling via equation Q does not exactly preserve the shape of 
the J ^ (J) tail, which corresponds to points very near to a single 
source; in this limit /(J) oc J^^^^ and all the curves should con- 
verge to the solid one. Fortunately, precisely because this tail de- 
scribes the immediate neighborhood of sources, it is not important 
for our argument; we are only concerned with the ionizing back- 
ground in the diffuse IGM, where transmission in the Lya forest 
appears. In practice, at the ionized fractions of interest most of the 
volume is filled by large bubbles (R > ro) for which the rescaling 
is quite accurate. 

More seriously, these curves assume that the galaxy density 
is Poisson-distributed around the universal mean value; in reality, 
galaxy clustering causes variations in the number counts above and 
beyond the Poi ssonian fluctuations in the halo number counts tha t 
we include (c.f. lWvithe & Loehll2006l ; lMesinger & Diikstrj|2008h . 
These have fractional amplitude ~ bia"(ro) in the fully-ionized 
limit, where bi is the mean bias of ionizing sources and cr^(ro) 
is the variance in the dark matter field smoothed on a scale ro ; for 
z = G and ro = 20 Mpc, we have bi(T(ro) ~ 0.25. To obtain /(J) 
throughout the universe, the solid curve in Fig.[T]must be convolved 
with the underlying distribution of N de termined by clustering, 
which broadens it substantially jMesinger & Diikstrall2008l) . For- 
tunately, this only strengthens our conclusions (see the discussion 
in so we defer a detailed model of it to future work. 

Finally, we also assume a spatially constant ro throughout the 
Universe. This is certainly an oversimplification, because ro must 
itself depend on J. However, the largest variations in J are across 
bubbles with R <^ ro, for which the attenuation length does not 
strongly affect our calculations. When ro is important, the varia- 
tions in J are relatively modest anyway, so our assumption seems 
a reasonable one. 



Figure 2. Distribution of J during reionization relative to its mean 
value in a fully-ionized IGM, (>.')ioni2e(j. The solid curves take Xi = 
0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 0.95, from left to right. The dot-dashed 
curve is for a fully-ionized universe. The dotted curve shows f(J) when 
Xi = 0.4 if we ignore fluctuations inside the discrete bubbles. All assume 
ro = 20 Mpc and z = 6.5. 



3 THE IONIZING BACKGROUND DURING 
REIONIZATION 

To obtain /(J) during reionization, we must convolve the dis- 
tributions for single bubbles in Figure [T] with the distribution of 
ionized bubble s izes, nf,(R). We use the analytic formulation of 
iFurlanetto et al.l ( l2004h to compute the latter; note that it is deter- 
mined primarily by t he ionized fraction and has only a weak de- 
pendence on redshift (Furlanetto et al. 20061; iMcOuiim et alj|2007l ; 
Mesinger & Furlanetto 2008), so our results are fairly generic. We 
also use this formalism to co mpute the mass functio n of dark mat- 
ter halos inside each bubble dFurlanetto et alj[2004) ; note that this 
procedure does properly capture clustering on scales R < ro, be- 
cause it is precisely this galaxy clustering which determines nb{R). 
In other words, overdense regions with more galaxies are already 
par t of larger ionized bubb les, with a correspondingly more galax- 
ies dFurlanetto et al.l2004h . 

Figure |2] shows /(J) at z = 6.5 for a variety of ionized frac- 
tions [normalized so that J dj f{J) — Xi, the mean ionized frac- 
tion]. In all cases we assume that r o = 20 Mpc, consistent with the 
extrapolations of iLidz et al.l ( |2007|) from models of the Lya forest 
at lower redshifts. The dot-dashed curve shows /(J) once overlap 
is complete. Again, the most important result is the relatively slow 
evolution of /(J) throughout the last half of reionization; (J) in- 
creases primarily because the tail toward small J shrinks as more 
and more of the Universe is incorporated into large bubbles. 

Because the ionized bubbles reach such large sizes early in 
reionization (typically ~ 4.5 Mpc when Xi — 0.5), the radiation 
background incident on a typical ionized patch increases by only 
a factor ~ 5 throughout the last half of reionization: although the 
size of each ionized region continues to increase until i? — > oo, the 
fixed attenuation length prevents a corresponding rapid increase in 
(J). Moreover, regardless of Xi, some regions have J ^ (J") ionized 
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Figure 3. Effective optical depth (a) and mean ionizing baclcground (h) as 
a function of Xi at a constant redshift. The filled squares, hexagons, and 
triangles assume ro = 10, 20, and 30 Mpc, respectively, at 2 = 6.5. 
The open hexagons show the mean ionizing background inside ionized re- 
gions, (r_i2) /xi, for ro = 20 Mpc at z = 6.5. The crosses use f(J) 
taken from the semi-numeric simulations of iMesinger & Diikstral feOOSh 
(see text). 



even though they are not particularly near any individual source. 
This broad range occurs because rn, {R) always contains some large 
bubbles, even early in reionization, and is not due to random Pois- 
son fluctuations in the galaxy counts inside each bubble. To see 
this, the dotted curve shows what happens if we ignore fluctuations 
within each bubble at Xi — 0.4 (i.e., we set J equal to its mean 
value for that R); it traces the full curve closely except in the high- 
J tail that is due to proximity to a single sourceQ 



4 OBSERVABLE IMPLICATIONS 

We now examine how the slow evolution of /(J) affects the Lya 
forest. The filled hexagons in Figure |3j> show the mean ionization 

rate (F) for the models in Figure l2Fl the filled squares and triangles sonable TbB « 8 at z = 6.5 (c.f. 



Figure 4. Same as Fig.[3}i, except for T^ff , 



rate of evolution decreases with rp , because bubbles reach the "sat- 
uration radius" R — ro quicker jFurlanetto & OhiboOSi) ; beyond 
that point, (J) can only increase by a factor (1 — l/e)~^ so long 
as ro remains fixed. 

The most easily observed property of quasar Lya forest spec- 
tra is the mean transmission T = exp(— Tcft); we show the ef- 
fective optical depth in Figure [3ji for the same set of models. To 
compute Tcff, we convolve the (volume-averaged) IGM density dis- 
tribution of Miralda-Escude et al. (2000) with /(J). We assume 
that the local density is uncorrelated with J; while this is certainly 
not correct in detail, it is a reasonable first step because F is typ- 
ically dominated by distant galaxies (Olber's paradox). Note that 
the absolute value of Toff is somewhat sensitive to the high- J tail; 
we truncate all of our calculations at 10 (J) j^^^ (this essentially 
prevents the low-density voids that provide transmission from be- 
ing extremely close to galaxies). For concreteness, we assume that 
(F_i2) = (F) /(lO"^^ s"^) = 0.05 when Xi = 1, consistent with 
limits from Lyg forest measurements in q uasar spectra at z ~ 6 
jFan et alj2002ll200^ : lBolton & HaehneltllT 



Fan et al 



2007,). This yields a rea- 



2001 



l2006l) . 



show the same sequence of ionized fractions, but with ro = 10 and 
30 Mpc, respectively (approximately spanning the range of values 
expec ted f rom extrapolating Lya forest measurements at lower red- 
shifts; iLidz et al. 2007). Note that these values are averaged across 
the entire Universe; the mean values within ionized regions are 
shown for the ro = 20 Mpc model by the open hexagons. Again, 
it is obvious that (F) evolves only slowly during reionization. The 



^ Imposing a duty cycle on star formation (and emission of UV pho- 
tons), or suppressing star formation in small halos, also has no substantial 
effect. Even if ~ 0.01, significantly smaller than one would expect from 
the dynamical times of galaxies, the shape of /(J) changes only slightly, 
except in the high- J tail (which increases in amplitude, because each galaxy 
dominates a larger region around itself). 

^ Here we assume F oc J, neglecting the variations in the mean free path 
for higher-energy photons because they do not contribute strongly to the 
total ionization rate. 



Figure |3] shows that there need not be a discontinuity in TcB 
near the end of reionization: in all of our models, the transmis- 
sion evolves smoothly and (relatively) slowly as a function of ion- 
ized fraction, especially when ro is small (so that bubbles enter the 
saturated regime earlier). However, the observational relevance of 
these results is questionable, because such large Toff are probably 
unobservable. Fortunately, the higher Lyman-series transitions are 
easier to observe; for example, Ly/3 has an optical depth 6.24 times 
smaller than Lya (although the ratio Tefi,/3 /r^ff i s closer to unity be- 
cause of the convolution with the density fi eld; ISongaila & Cowij 
l2002l ; ISongailj2004l ; IOh & Furlanettcl2005l) . FigureHshows r^s^p 
for the same models as Figure [3}i. This shows the same qualitative 
trends as for Lya, although the evolution is somewhat steeper, and 
again illustrates that no strong break is necessary at overlap. 

As described above, the most important shortcoming of our 
model is that it ignores galaxy clustering when R > ro; cluster- 
ing substantially broadens /(J) (but does not affect (J) much). To 



© 0000 RAS, MNRAS 000, 000-000 



The ionizing background during overlap 5 



check the importance of this effect, we have drawn ,f(J) from the 
z — 7 semi-numeric simulations of 'Mesinger &. Dijkstra ('20 St) 
(in tu rn based on the reionization model of Mesinger & Furlanettal 
l2007h . These simulations include attenuation (with ro = 20 Mpc), 
clustering, and more accurate bubble sizes. In principle, we can 
also use them to compute correlations between the flux field and 
the underlying density field; however, the finite box size (100 Mpc) 
means that rare voids (which account for the bulk of the transmis- 
sion) may be missed and the quasi-linear treatment may not capture 
the full density distribution. We therefore dra w f(J) from the semi- 
nume ric simulations and convolve it with the lMiralda-Escude et al.l 
1 I2OOO ) density distribution to compute Tcff, shown by the crosses 
in Figures [3] and I4FI These show even slower evolution through- 
out reionization, partly because the simulated ionized regions are 
somewhat larger than the analytic model predicts and partly due to 
the inclusion of clustering when R > ro. Thus the semi-numeric 
approach confirms our conclusion that overlap need not be accom- 
panied by a rapid increase in (T) , unless ro also increases. 

To this point, we have fixed the mean free path and redshift 
to isolate how Xi affects /(J). With the filled hexagons in Fig- 
ure [5] we specialize to a particular reionization scenario where 
Xi is proportional to the fraction of gas in galaxies, and where 
we set Xi = 1 at 2; < 6.5. We further assume that ro — 
15[(1 + z)/7.5]~'^ Mpc, consistent with extrapolations from the 
lower redshift forest jLidz et alj2007t) . Again, we normalize so that 
(r_i2) = 0.05 at 2; = 6.5. Even with this evolution in ro, the mean 
transmission evolves smoothly across this entire range, and there is 
only a modest break in (T) at z — 6.5. 

In contrast, the open squares assume Xi — 1 at all redshifts, 
with an identical ro{z) and a constant comoving emissivity. The 
latter does give more transmission at all redshifts, but both models 
evolve quite smoothly - there is only a slight break in the reion- 
ization model's slope at overlap compared to the fully-ionized one. 
Thus overlap does not, in general, cause an obvious feature in the 
transmission, even for the higher-series transitions. In fact, slow 
evolution is guaranteed to occur unless the mean free path itself 
evolves more rapidly during reionization. 



5 DISCUSSION 

We have shown that, because attenuation due to LLSs must become 
important before overlap, the late stages of reionization need not be 
accompanied by a rapid increase in (F) (or Tcff). Thus, the conven- 
tional wisdom that a sharp increase in the ionizing background is a 
robust indicator of overlap must be modified: namely, such a feature 
indicates only that the attenuation length ro is evolving rapidly. 

Of course, it may be that ro does evolve quickly at the end 
of reionizat ion (and some theoretical calculations sugg est that this 
does occur; IWvithe et alj|2008l; Ichoudhurv etai]|200a) . For exam- 
ple, (r) must (at least to so me degree) control the ionization struc - 
ture of LLSs and hence ro dMiralda-Escudi |2005| ; ISchavell200^ . 
As (r) increases, LLSs will shrink, increasing ro and hence (F), 
creating a positive feedback loop. Thus even a slow initial increase 



We actually take the 2 = 7 f{J) from the simulations and apply them 
to the z = 6.5 density field. We then choose (r_i2} = 0.032 at = 1 
so that Teff is identical to that in the analytic model. Note that, by widen- 
ing /(J), clustering decreases the required mean ionizing background for 
a given mean transmission level. Interestingly, it also changes the relative 
ratio of Lyo and Ly/3 transmission. We have verified that the slow evolution 
persists even if we use the semi-numeric density field. 
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Figure 5. Evolution of the effective optical depth in the Ly/3 transition (a) 
and the mean ionizing background (b) for a scenario in which reionization 
ends at z < 6.5 (marked by a vertical dotted line) and ro oc (l + z)"^ 
(filled hexagons). Above 2 = 6.5, adjacent hexagons have Axi = 0.1 
between them (with Xi = 0.5 at z = 7.8). The open squai'es show an al- 
ternate model with a constant comoving emissivity and Xi = 1 throughout 
this range. 



in the emissivity and/or mean free path at the end of reionization 
may spiral into a relatively rapid change in (F). Unfortunately, our 
understanding of LLSs is not yet advanced enough to determine the 
effectiveness of such a feedback loop; for now, we can only say that 
effective feedback requires that th e LLSs have quite stee p density 
profiles (see, e.g., the Appendix to lFurlanetto & Ohll2005l) . 

More importantly, there is no obvious reason that the mean 
free path can only change during overlap: so long as (i) the emis- 
sivity continues to increase and (ii) (F) controls the properties of 
LLSs, the feedback loop will continue - and with the rapid increase 
in the collapsed fraction at « > 6, such a scenario see ms entirely 
plaus i ble. T hus, even if (F) is evolving strongly at 2 ~ 6 jFan et all 
I2002l l200g), we do not have a robust indicator of overlap. 

Indeed, it may be better to regard the final stage of reionization 
as the disappearance of LLSs after overlap. This "post-overlap" 
phase, which consumes just a few percent of the neutral gas, is typ- 
ically thought to be indistinguishable from the cosmic-web domi- 
nated Universe at later times. We have shown instead that (F) may 
evolve slowly during overlap but rapidly afterwards, making quasar 
absorption spectra a useful probe of this tail end of reionization 
(which matches smoothly onto the post-reionization Universe) but 
not necessarily of overlap itself. This contrasts with the conven- 
tional picture in which (F) evolves rapidly only during overlap; that 
intuition came from reionization simulations that were too small to 
include the full inhomogeneity o f reionization and so exaggerated 
the importance of overlap (e.g., Gnedin '2000") . The recent semi- 
numeric calculations of Choudhury et al. (2008), which included 
an approximate prescription for self-shielding gas, present a sim- 
ilar picture to ours, in which these LLSs are burned off over the 
(rather extended) final stages of reionization that follow overlap. 

Beyond an increase in F, there are other reasons that ro may 
evolve throughout (and after) reionization. For example, photo- 
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heating can evaporate loosely-bound structures, substantially mod- 
ifying the IGM gas distribution by eliminating dense systems 
dHaiman et al. 2001; Shapiro et al. 2004; Pawlik et al. 2008). This 
will increase ro as well; however, the process should happen grad- 
ually throughout reionization as more and more of the volume 
is illuminated (and so will be particularly slow if reionization is 
extended). Evaporation also occurs over the sound-crossing time, 
which is relatively long for the moderate overd ensities of most in- 
terest at these high redshifts jPawlik et alj200^ . so that LLSs prob- 
ably continue to evolve past the "end" of reionization. 

Given the challenges we have raised to the conventional inter- 
pretation, it is worth asking two further questions about the data. 
First, are there other effects that can cause a rapid increase in TcS 
without a significant change in F? One possibility is the density 
distribution itself: at least according to current models, it is ex- 
tremely steep in the low-density tail that is responsible for Lya 
forest transmission, so it is possible that a small increase in the 
ionizin g background renders a relatively large fraction of the IGM 
visible loh & Furlanettol2005l) . This must await more detailed cal- 
culations of the evolving density field at the end of reionization. 

Second, how certain can we be that overlap has actually oc- 
curred by 2 ~ 6? Is it possible that, although most of the IGM is 
highly ionized before that point, a small f raction far from ionizing 
sources could still be completely neutral l lLidz et aT]|2007l) ? In the 
conventional picture, in which (F) evolves rapidly at overlap, such 
a conclusion can be easily dismissed. However, our results suggest 
that overlap itself may be buried inside the smoothly evolving trans- 
mission at 2 < 6 - or it may have occurred at much higher redshifts. 

Finally, we expect our conclusions to hold with even more 
force during helium reionization: the clumpy IGM and enhanced 
recombi nation rate during tha t era makes attenuation more im- 
porta nt jFurlanetto & Ohll2008l ; iMcOuinn et alJl2008l : lBolton et al.l 
l2008l) and high -energy photons can more easily create a nearly- 
uniform ionizing background, so the transition from bubble- 
domination to web-domination will be even smoother. 

In summary, with our present knowledge of high-2 LLSs there 
is no compelling reason to associate evolution in ro or Lya for- 
est transmission exclusively with overlap. If inde ed the quasar dat a 
shows a more rapid evolution in Teg at 2 > 6 jFan et al]|2006l) . 
this may indicate that we are seeing a rapid increase in ro that 
followed or even preceded ove rlap (note, however, that this mea- 
surement is itself controv ersial; ISongaila & Cowij|2002l ; ISongailal 
l2004lBeckeretal.ll2007l) . More detailed studies of the coupling 
between the ionizing background and the dense, neutral blobs that 
trap ionizing photons are needed before the next step - associating 
such a jump with overlap itself - can be taken securely. 
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